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Abstract: We extend the Sakai-Sugimoto holographic model of QCD (HQCD) by including 

the scalar bi-fundamental "tachyon" field in the 8-brane-anti-8-brane probe theory. We show 

that this field is responsible both for the spontaneous breaking of the chiral symmetry, and for 

the generation of (current algebra) quark masses, from the point of view of the bulk theory. 

As a by-product we show how this leads to the Gell-Mann-Oakes-Renner relation for the pion 

H ' mass. 
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1. Introduction 

The closest model so far to a holographic description of large N c QCD is the Sakai-Sugimoto 
model 0. The U(N C ) gauge sector is described, at low energy, by the near-horizon limit of 
N c D4-branes wrapped on a Scherk-Schwarz circle @. This description is valid at energies 
well below the Kaluza-Klein scale of the circle. The quark sector is incorporated by includ- 
ing Nf D8-branes and Nf anti-D8-branes transverse to the circle. The strings that stretch 
between the original D4-branes and the D8-branes (anti-D8-branes) describe right-handed 
(left-handed) chiral fermions, which transform in the fundamental representation of both the 
U(N C ) color group and the U(Nf)n (U(Nf)i,) flavor group. It is assumed that Nf <C N c , 
so that the D8-branes can be treated as probes in the D4-brane background, and one can 
ignore their backreaction. In QCD this corresponds to the quenched approximation. The 
Sakai-Sugimoto model shares many features with other holographic models of gauge the- 
ory with matter, however the novel feature of this model is the geometrical realization of 
spontaneous flavor chiral symmetry breaking. The 8-branes and anti-8-branes are separated 
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Figure 1: The 8-brane-anti-8-brane configuration in the compact non-extremal 4-brane background, 
(a) The anti-podal case, (b) The general case. 

along the circle asymptotically in the radial coordinate, but are connected at some minimal 
radial position (figure [l]). The former corresponds in QCD to the UV flavor symmetry be- 
ing U(Nf)n x U(Nf)i, and the latter corresponds to the IR symmetry being the diagonal 
subgroup U(Nf)v- 

In spite of its success, there are several open questions about the model, some of which 
are related to very basic notions of gauge dynamics. The first is the incorporation of a QCD, 
or "current algebra" quark mass. The quarks are massless in this model since the 8-branes 
necessarily intersect the 4-branes. This is also manifested in the fact that the modes identified 
with the pions are massless. It is well known that pions obey the Gell-Mann-Oakes-Renner 
relation 101 



M|=« + «). (1.1) 

J IT 



Therefore in a non-trivial quark condensate massless pions imply massless quarks. In the 
original construction of Q the 8-branes and anti-8-branes were located at antipodal points on 
the circle, and they connected at the minimal radial position of the background ukk- This was 
extended in 0] , || , |g] to a family of configurations parameterized by the asymptotic separation 
L, or equivalently by the minimal radial position of the 8-branes uq (see figure 0(b)). In these 
configurations there is a natural mass parameter associated with the "length" of a string 
stretched from the minimal radial position of the background ukk to the minimal radial 
position of the 8-branes uq. However, it is easy to check that the pions remain massless, so 
this parameter cannot be identified with the current algebra mass. It is instead related to 
the "constitutent quark mass". Indeed in || for the model of [|J, and in Q for analogous 
non-critical models, it was found that the masses of the vector mesons are linearly related 
to this "length of the string" mass parameter. Moreover, a model of the decay process of 
spinning stringy mesons Q supports the interpretation of the "length" as the constituent 
mass. 

A related question has to do with the quark condensate itself: how does one compute 
it in the holographic description? 1 The answers to both questions are related to each other. 

1 In |8| both the current algebra mass of the quarks and the condensate can be read from the profile of the 
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The quark mass term in QCD is 



m q qq = m q (q R q L + q[q R ) , (1.2) 

where m q is an Nf x Nf matrix, g# is a fundamental of U(Nf)n, and qi is a fundamental 
of U(Nf)i. Both the mass and the quark bi-linear should therefore be identified with a bi- 
fundamental field in the bulk. 2 In this model the required bi-fundamental field comes from 
the D8-D8 strings. 3 The dual operator is therefore non-local in the coordinate transverse to 
the 8-branes. According to the usual holographic dictionary the normalizable mode of this 
field should correspond to the expectation value of the quark bi-linear, i.e. to the quark 
condensate, and the non-normalizable mode should correspond to the quark mass. 

The bi-fundamental field was not included in the analysis of Q, since, as was argued 
there, it is very massive. The 8-brane-anti-8-brane separation was assumed to be much 
greater than the string length, which in flat space would make this field massive. However 
the proper distance between the 8-branes and anti-8-branes in the curved background of this 
model depends on the radial coordinate u, and decreases as u decreases. The mass of the 
bi-fundamental field therefore depends on u as well. For the U-shaped configuration found in 
[0 this field remains massive for all u. While the proper distance decreases as u decreases, 
the relative angle between the 8-branes and anti-8-branes increases, so that it never becomes 
tachyonic. The situation changes in the case of the non-compact background considered in 
Q. In that case ukk = 0, so there are two possible configurations: a connected U-shaped 8- 
brane similar to the one of Q, and a disconnected parallel 8-brane-anti-8-brane configuration. 
In the parallel configuration the proper distance between the 8-branes and anti-8-branes goes 
to zero at u = 0, and therefore the bi-fundamental field becomes tachyonic in a finite range of u 
near the origin. This represents a (radially) localized tachyonic instability, and one expects the 
true vacuum to be the U-shaped configuration 4 . The condensation of this localized tachyon 
can be seen as a Higgs-like effect which breaks the chiral symmetry U{Nf)n x U{Nf)i to the 
diagonal group U(Nf)y. 

In either case, the remaining non-tachyonic mode of the bi-fundamental field, which we 
will continue to call the "tachyon" T, is crucial for describing the quark mass and condensate. 

In this paper we incorporate the tachyon into the 8-brane action using a proposal of 



Garousi for the brane-antibrane effective action [14], which extends Sen's original proposal 



for the non-BPS D-branes [15|. We show that the coupled equations of motion for the tachyon 



T and the 8-brane-anti-8-brane separation L admit a solution which describes a U-shaped 



flavor branes, and the GOR relation of ([Tl]) is obeyed. However that model suffers from the drawback that it 
does not incorporate chiral flavor symmetry. 

2 ln U the quark mass and condensate were identified with the scalar field corresponding to the 8-brane- 
anti-8-brane separation. We believe this is incorrect, since this field transforms in the adjoint, rather than the 
bi-fundamental, representation of the U(Nf)s. 

3 The same field appears in the holographic description of the resolution of the (7(1) a puzzle E3 |. 

4 A similar effects occur in the "hairpin brane" of [Q, and in the meta-stable supersymmetry breaking 
brane configurations of mM. 
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configuration. In our solution the tachyon has a non-trivial profile, which for large u is a 
linear combination of a normalizable mode and a non-normalizable mode. We relate the 
coefficient of the former to the quark condensate (qq), and the coefficient of the latter to 
the quark mass m q . We also show that the pions, which are part of the meson spectrum, 
acquire a mass that satisfies the GOR relation (1.1). For m q = our solution describes the 
same configuration as [ll]] , but it also includes the effect of the (normalizable mode of the) 
massive bi-fundamental field. At large u the solutions are the same, but the precise shape of 
the 8-brane at finite u changes. 



For simplicity we will consider the non-compact case dual to the NJL model [|ll| , namely 
the near- horizon background of N c extremal D4-branes. The metric in this case does not 
contain the "thermal factor" f(u) = 1 — v? KK /v?. The behavior near the boundary at u — ► oo 
will be similar to the compact case since /(oo) = 1, and therefore our results for the quark 
mass and condensate, which are determined by the behavior of T near the boundary, will be 
the same. This is also reasonable from the field theory point of view. While the gauge sector 
of this model is very different from QCD, and Kaluza-Klein states do not decouple, the flavor 
sector, which is where chiral symmetry breaking and quark masses are seen, is the same. We 
will also deal only with the one flavor case Nf = 1, for which the 8-brane theory is Abelian. 
Note that in this case the would-be broken symmetry is the anomalous U(1)a- At large N c , 
however, the anomaly, and with it the mass of the would-be Goldstone boson r/, is suppressed 
[|l^ , |l§| l (see however [10] for a discussion of how it is suppressed in this model). The GOR 



relation (1.1) therefore holds also for the rf . 
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A holographic dual description of the chiral condensate and quark mass in QCD has been 



discussed previously in the context of the "bottom-up" AdS/QCD model |2(], 21], which is es 



sentially a five-dimensional U(Nf) x U(Nf) Yang-Mills theory in AdS§ with a bi-fundamental 



tachyonic scalar field. This was later generalized to a tachyonic DBI + CS theory in [22]. 
The outcome of the present paper is a holographic picture where 

• The spontaneous breaking of flavor chiral symmetry emerges from a "Higgs mechanism" 
with an order parameter which is the expectation value of the bi-fundamental tachyon 
field. 

• The current algebra mass of the quarks is associated with a non-normalizable mode of 
the tachyon. The quark anti-quark condensate can be identified with a normalizable 
mode of the tachyon. 

• The pions of the model obey the GOR relation. 

The paper is organized as follows. In section 2 we review the proposal for the Dp-Dp 
effective action, and apply it to the D8-D8 system in the near-horizon extremal 4-brane 



5 A different approach to the pion mass in this model was discussed in Jig ]. 

6 There is an alternative large N c extension of one-flavor QCD, in which the fermions transform in the 
anti-symmetric representation of the gauge group j^]. In that model the anomaly is not suppressed, and the 
GOR relation is not expected to hold for the rf . 
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background. In section 3 we study the asymptotic forms of the solutions for T{u) and L(u), 
both at large u and near the point u = uq where the branes and anti-branes connect. We 
extract the quark mass and condensate from the behavior of T at large u. In section 4 
we present numerical solutions which interpolate between the two asymptotic solutions, and 
compare with the solution without the tachyon of flllfl . In section 5 we begin to analyze the 
meson spectrum in the tachyon background. This includes both the fluctuations of the scalar 
fields T and L, as well as the worldvolume gauge fields on the 8-branes and anti-8-branes. In 
particular we show that the mass of the pions satisifies the GOR relation. 



2. The D8-D8 theory 



A proposal for the effective action of a parallel p-brane-anti-p-brane system in curved space- 
time was given by Garousi in [14|. Denoting by and the adjoint (position) scalar 
fields and gauge fields on the branes (n = 1) and anti-branes (n = 2), and by T the complex 
bi-fundamental scalar field, the action is given by 



S = -T p f <F +1 a e-^ x( ^V(T)Jl + 

71=1,2 



*(* (n) )T/m«/i + |T|2|L| \ /-det(g(")+TH), (2.1) 



where L = — X^ is the brane-antibrane separation, and 



rr-in) 
1 ab 



p(n) 



G, 



ab 



G nn UVG 



+ 2na'F, 



1 + 



|T| 2 |L| : 



I rp{n) 



ab 



ira' (D a T(D b T)* + D b T{D a Tf) 



+ I (G ai + d a X^Gji) V (T(D b T)* - T*D b T) 



+ % - (T(D a T)* - T*D a T) V (Ga - G tj d b X^ 



(2.2) 



(2.3) 



We use a, b for the worldvolume directions, and i,j for the transvese directions. The covariant 
derivative of the bi-fundamental scalar is given by D a T = d a T - - A^)T, and V(T) 

is the scalar field potential. 

This action was obtained by generalizing Sen's action for a non-BPS 9-brane in Type IIA 
string theory |19] as follows. First, the tachyon kinetic term is added under the square root 
p5|]. Second, the action is extended to two unstable 9-branes by the familiar symmetric trace 
prescription for the non-Abelian DBI action. Third, the action is transformed to a 9-brane- 
anti- 9-brane action in Type IIB string theory by projecting with (—1) Fl . Finally, the general 
p-brane-anti-p-brane action is obtained by T-duality. To separate the branes and antibranes 
we turn on a Wilson line in the 9-brane-anti-9-brane model, which fixes the dependence on 
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Figure 2: The inverse cosh Tachyon potential. 



L. As a check, note that for L = this reduces to Sen's action for a coinciding brane and 
antibrane [p4[] . 

The tachyon potential for the brane-antibrane pair is not known precisely even in flat 
space. Boundary superstring field theory gives a potential |2j| 

V (T) = e~ T2 ' A . (2.4) 



An alternative proposal for the potential is [p6| 27, 28, 14] 

^-ssjW <2 - 5) 

This reproduces, for example, the S-brane solution using the tachyon effective theory |pq |. In 
both proposals (and there may be others) the true vacuum is at T — > dboo, but the details 
are different. We will work with the inverse cosh potential (^^) shown in figure ^. 

Let us apply this proposal to the D8-D8 system in the non-compact extremal D4-brane 
background. The background is defined by 

3 _ 3 

ds 2 = 2 [r }flv dx»dx v + (dx A ) 2 ) + (^j 2 {{dUf + (Ufdnj) , (2.6) 



3 



and 

«*-ft^J ■ (2-7) 

where n,v = 0,1,2,3, and < U < oo. We assume that the 8-brane and anti-8-brane 
are positioned symmetrically at = L/2 and = —L/2, respectively, and that the 
configuration depends only on the radial coordinate U. It will also be convenient to work in 
the unitary gauge in which T is real. Suppressing the gauge fields for now, the action for the 
D8-D8 pair in this background becomes 

S[T,L] = - 2M j d 4 xduV(T)u 4 ^V[T, L] , (2.8) 
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where 

V[T,L] = 1 + ^(L'(u)f + ^(T'(u)f + -J-^( L (u))\T(u)f , (2.9) 

and where we have defined u = U/R and J\f = ^Q^R 5 jg s . 

Note that the proper distance between the 8-brane and anti-8-brane is 

^proper — 

u 3/4 L, (2.10) 

so even if we keep the brane and antibrane well separated in coordinate distance, the proper 
distance will decrease below the string scale for small enough u, and the field T will be 
tachyonic in that region. One can see this directly by expanding the action for small T, 
which gives (after properly normalizing to get a canonical kinetic term) 

m|w = __l 7 + (^W)!. (211) 

We recognize the first term as the zero-point energy of the open superstring in the NS sector 
in flat space, and the second term as the contribution of the proper length of the open string. 
This result is most likely not precise. First, the flat space result for the zero-point energy 
probably changes in this background. We do not know how to compute it, since this is an 
RR background. Second, the straight string stretched between the 8-brane and anti-8-brane 
is not the minimal length (and mass) string, and it actually prefers to curve down in u [11]. 
However we believe that the qualitative result is still correct, namely that m\ < below 
some critical u. In other words the field T has a localized tachyonic mode in a small region 
near u = 0. This is similar to the tachyon which appears at the intersection of branes which 
meet at a small angle. We therefore expect the ground state to correspond to the connected 
8-brane configuration. 

2.1 The compact case 

In the Sakai-Sugimoto model x± is compact and the near-horizon metric is given by 

3 3 
ds 2 = (I) ' (v,udxV + f(U)(dx A f) + 5 (£QjL + V 2 dn 2^ (212) 



where 



m) = l-%t,U KK = l*. (2.13) 



Strictly speaking, the 8-brane and anti-8-brane cannot be treated as separate entities in this 
background, since the X4 circle shrinks to zero size at U = Ukk- In this case the brane and 
anti-brane are necessarily connected, and the D8-D8 action should be viewed as a large-u 
effective theory for the worldvolume fields on the two sides of the 8-brane, together with the 
massive scalar field coming from the open string stretched between the two sides. At large U 
the compact background is essentially identical to the non-compact one, so the results related 
to the mass and condensate will be the same. The precise profiles of the fields T(u) and L(u) 
at finite u will be different. 
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3. Asymptotic solutions 



The equations of motion that follow from Q2.8|) are given by 



d_ 

du 

d_ 

du 



V(T) u 4 



V(T) ui , 



L' 



V{T) 



V(T) 



V 



u?LT 



dT 



(3.1) 
(3.2) 



where T> was defined in (|2.9|) , V{T) is the inverse cosh potential, and we have set 2-7rc/ = 1. 

The tachyon equation (|3.2[ ) has a trivial solution T = 0. In this case the solution to the 
L equation Q3,l| ) is L(u) = L^, corresponding to the parallel D8-D8 configuration. 7 This 
configuration is unstable due to the localized tachyon mode near u = 0. The stable solution 
must involve a non-trivial tachyon condensate T(u), which, as we shall see below, corresponds 
to a single U-shaped configuration. 

We expect the 8-brane and anti-8-brane to connect roughly at the radial position below 
which the bi-fundamental field is tachyonic. Let us first expand the fields near this point: 



L{u) = (u - u ) p [l + h(u-u )-\ ] 

T(u) = (u - u ) q [t + t\{u — uq) H ] , 

where we assume that lo,to > 0. To leading order, the L equation fl3.1|) gives 



-2 and to 



2R 2 



K 3/2 



(3.3) 
(3.4) 

(3.5) 



This implies, in particular, that p > 0. In addition, the absence of sources at u = uq implies 
that p < 1, so that L'{uq) — > oo and the configuration is smooth. The T equation (|3.2| ) is 
then also satisfied to leading order. The leading behavior near uq is then 



T(u) ~ (u - uq)' 



L{u) 



(U 



< p < 1 . 



(3.6) 



This is in accord with the interpretation of the non-trivial solution as the chiral-symmetry- 
breaking U-shaped configuration. The brane-antibrane separation vanishes at u = uq, and 
the tachyon diverges, i.e. goes to its true vacuum in the potential (f2~~5|). 

To compute the gauge theory quantities, in this case the quark mass and condensate, we 
should look at the behavior of the solution at large u. This corresponds to the UV limit of 



7 Equation (3.1) reduces in this case to the same equation one gets from the single 8-brane action without 
the tachyon Jll|. There are two solutions in that case corresponding to a straight 8-brane and a U-shaped 
8-brane. However the action in our case is doubled since it includes both an 8-brane and an anti-8-brane. 
Consequently there are four possible solutions with T = 0, corresponding to either brane or antibrane being 
straight or U-shaped. We are interested only in solutions with two asymptotic boundaries, one for the 8-brane 



and one for the anti-8-brane (or equivalently we require X^\u) and Xf'{u) to be single- valued) . With T = 
that leaves only the straight and parallel D8-D8 solution. 



- 8 - 



the gauge theory. Strictly speaking, the UV limit is not well-defined in this model, since it is 
really a five-dimensional gauge theory. We will therefore always be considering a UV cutoff 
Uoq. In this regime the field T is very massive, so we can consider small fluctuations away 
from the trivial solution: 



L(u) = Loo + L(u) 
T(u) = + f{u) , 

where L <C and T< 1. In this approximation the action is quadratic 



S oc / d A xdu 



u^ 2 + 



1 



8R 2 



11/2/ f,,2 



_l_ ( f/ )2+ ^o (f) 2 



The asymptotic solutions for u>l are given by 

L(u) C L u- 9 ' 2 



f{u) 



u- 2 (C T e 



-RLooU 



+ C' T e 



OC U 



(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 



These solutions are only valid in a regime of u for which the fluctuations are small. We 
therefore have to assume that C' T < e -R L cx>v-oo^ w hereas Ct and Cj_, can be taken to be C(l) 
in the cutoff. 

3.1 Quark mass and condensate 

The growing and decaying exponentials correspond to the non-normalizable and normalizable 
solutions for T, respectively. We would therefore like to identify the coefficients C' T and Ct 
with the quark mass m q and quark condensate (qq), respectively. Let us verify this explicitly. 
In QCD (at zero temperature) the quark condensate is given by the variation of the energy 
density with respect to the quark mass 

SSqcd 



(w) 



5m n 



Let us assume that m q is given by the (dimensionless) parameter C' T 

AC T , 



(3.12) 



(3.13) 



where A is some fixed mass scale. To evaluate (3.12) in the holographic dual we must vary the 
asymptotic (Euclidean) 8-brane action ( |3.9j ) with respect to the parameter C' T of the solution. 
The general variation is 



5S 



du 



< )C 5 T{u) + -^-ST'(u) + J£-5L(u) + -J^-rSL'(u) 



_6T(u)" v ' ' 5T'(u 
Using the equations of motion this reduces to 

5C 



5L(u) 



5L'{u) 



5S 



5T'{u) 



5T(u) 



+ 



_5C_ 

5L'(u 



■8L(u) 



(3.14) 



(3.15) 
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Focussing on the variation with respect to the tachyon we find 

2M u 5 / 2 V(T)T'(u) 



5S 



R 2 



V 



-8T(u) 



(3.16) 



Only the upper limit contributes, since although both T and T' diverge in the lower limit, 
the potential V(T) ~ exp(— y^T) — ► much faster. Using the large u asymptotic form of 
the solution ( 3.11] ) we find for a variation with respect to C' T : 



5S 

Jo? 



2M 



Loo {Ct — C' T e 



) 



(3.17) 



where we have imposed the cutoff Uqq. Since C' T is identified with the quark mass, we find 
that the quark condensate is related to Ct as 



(qq) 



2ML 
AR 



(3.18) 



4. Numerical solutions 



The asymptotic solutions near u = uq and at large u must connect in the full solution to the 



equations of motion ( p.l| ) and (3.2). In this section we present a numerical analysis of these 
equations. For convenience we define the dimensionless quantities (recall that 2ira' = 1) 



y = 4R 4 u, f(y) = 
In terms of these the D8-D8 action (El 

N 



^L(u), g(y)^V2T( 



(4.1) 



becomes 



d xdy- 



y 



64i? 14 J cosh(y^g(y)) ' 

where 

v = y~ 3 + f{y? + y- 3/2 g'(y) 2 + y- 3/2 f(y) 2 g(yf 

and the equations of motion become 



d_ 

dy 

d_ 

dy 



y4V-l/2 f ,( y) 



cosh( v /f5(y)). 

y^V-^g'jy) 
cosh( A /f g(y)) 



cosh (\/fff(y)) 

yWT>-Wf(y)*g(y) 
cosh( v / f g(y)) 



tanh( 



g(y)) 



2 cosh( v ^< 7 (y)) 



y 4 vi 



(4.2) 

(4.3) 

(4.4) 
(4.5) 



The range for y £ [0, oo) will be approximated numerically by the range [0.01, 100]. 

The solution is fixed by imposing boundary conditions for f(y), g(y), and their derivatives, 
either at infinity (UV) or at y = 1 (IR), which corresponds roughly to uq. Let's look at UV 
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Figure 4: (a) The shape of the 8-brane, and (b) the profile of the tachyon from the IR. 



boundary conditions first. Guided by the UV asymptotic form of the solution (3.10) and 



( 3.11 ), we impose 

/(100) = 1, /'(100) = 10~ 10 , 5 (100) = 10~ 30 , s'(lOO) = -10~ 30 . (4.6) 

Figure ^(a) shows the resulting numerical solution for the shape of the 8-brane, and figure 
||(b) shows the tachyon profile for this solution. The tachyon increases as y decreases, and 
blows up, i.e. attains its vacuum value, where the 8-brane and anti-8-brane connect. 

Now let's look at IR boundary conditions. Guided by the IR asymptotics ( |3.6[ ) we impose 
the numerical boundary conditions 

/(I) = 0.001, /'(1) = 500, 5 (1) = 400, g'(l) = -16000. (4.7) 

The solution near the connection point at y = 1 is shown in figure The behavior is 
qualitatively the same as with the UV boundary conditions. When we look at larger y, 
however, the qualitative behavior changes (figure ||). This is due to the sensitivity of the 
numerical solution to the IR boundary conditions (|4.7| ). The exact IR boundary values of 
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Figure 5: Same as figure 3, but going to larger y. 




y 



Figure 6: The solid line describes f(y) in our tachyonic model. The dashed line describes the AHJK 
solution /ahjk(^)- 



f',g and g' are infinite. As we increase the numerical IR boundary values of f',g and g' the 
region where the behavior changes moves to larger and larger y. 

4.1 A comparison with the AHJK solution 

Let us compare this solution with the solution for the single 8-brane action without the 
tachyon found in |Q. The latter is just the solution of ( fl.4| ) with g(y) = 0: 



d_ 

dy 



2/ 4 /ahjk(^) 

Vv-' s + W(y) 



0. (4.8) 



Using the same UV boundary conditions we arrive at the numerical solution presented in 
figure ||, where we also present our solution for comparison. The two configurations are very 
close in this regime, which is understandable since the correction due to the very massive field 
T is small. 
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We would like to argue that our solution with the non-trivial tachyon profile is, in some 
sense, a better approximation to the exact string theory solution. As evidence for this we will 
show that the free energy of our solution is smaller than that of the solution of (llfl , which 
we will refer to as the AHJK solution. 8 The free energy is given by the Euclidean action of 
the solution, which for our solution is 



pr ,, * M f J vWy- 3 + f'iv) 2 + y- 3/2 g'(y) 2 + y-^ 2 f{y) 2 g{y? fAQ . 

The free energy for the AHJK solution is given by 

£ahjk [/ah jk (y )] = j ^V^r 3 + /ahjkO/) 2 • (4.10) 

By splitting the integration over y to a UV part and an IR part, we can estimate the free 
energy of our tachyonic solution and the AHJK solution. In the IR region, the free energy 
of the tachyonic solution £ IR is strongly suppressed by the exponentially vanishing factor 
cosh -1 (y / 7r/2g(y)), which comes from the tachyon potential. So £ IR is obviously much smaller 
than the IR part of the AHJK free energy £ahjk- On ^ e °ther hand, in the UV region, we 
shall compare £ and £ahjk numerically. From figure ^(b), and using the numerical solutions 
of /, g, /ahjk, we calculate 9 

64i? 14 uv uv 

\ c ~~ c AHJK J 



100 aWy- 3 + / /2 + 2/- 3 /V 2 + ir 3 / 2 /V f 100 , 4 r^TT — * 
a y r? jw \ / a yy yy +Jahjk 

20.7859 COSh {^J^g) 720.3596 V 

= -818.417 < 0, (4.11) 
that is to say, £ uv < <?ahjk- Combining the results in the IR and UV regions, we see that 

£ < £ahjk , (4.12) 
so our tachyonic solution appears to be more favorable than the AHJK solution. 



5. The meson spectrum 



We now turn to the analysis of the spectrum in the tachyon U-shaped background. The 
fluctuations of T, L and the gauge fields on the 8-brane and anti-8-brane correspond to various 
mesons, including scalars, pseudo-scalars, vectors and axial- vectors. We are interested mainly 
in the lowest pseudo-scalar modes (the pions). In particular, we would like to see how they 
acquire mass when the quarks are massive. For completeness we also set up the eigenvalue 
problems for the other mesons, but we leave the (numerical) analysis for future work. 

8 Recall also that the U-shaped solution of |n]] has a lower free energy than that of the parallel configuration. 
9 The lower bounds of the integration intervals depend on the capacity of our computer. 
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5.1 Scalar fields 

Let's start with the scalar fields L and T. We expand around the classical solution 

T(a?,u) = T{u) + t(x», u) (5.1) 
u) = L{u) + £(x^,u) , (5.2) 

where t(x^, u) and i{x^, u) are real scalar fields 10 . Expanding the 8-brane-anti-8-brane action 
to quadratic order gives 



S[t,l] =-M"J d 4 xdu 



Ji(^t) 2 + i 2 {d^y + Tipton + T 4 (0 + ^ + X 6 tt' 



where the coefficients are given by 



(5.3) 



V 4R 2 (1 + Jl 2 T 2 ) J ' ' V 4 4i? 2 (l + nlL 2 T 2 ) 

y(T) ^L'T' _y(T)/n§ u(T') 2 \ 

3 ' ^21P(1 + uS L 2 T 2 ) ' 4 ' ^ 2 ^j' 

V(T)uL 4 T 2 V(T)JL 2 dV(T) 2ul L 2 T d 2 F(T) 4 /- 

2-5 : = 3 1 7= 1 ™ 1= 1 ^ — u W, 

£>§ VV dT VP dT 2 

_ 2V{T)uL 2 TT dV{T) 2uiT' _ V(T)u 4 V(T)u 4 (L') 2 

6: ~ B^vl + dT R2^V> 7 : " 4R 2 VV ~ I6i?4pf ' 

_ V(T)u^T 2 V(T)uL 2 T A _ V(T)u^LL'T 2 V(T)uH'T 
M :— /= s > ^9 :— — 3 > Mo '■— ■ — 3 > 

5 



2V{T)uLT 2 T' V(T) U 2L 2 L'T dV{T) u A L 

Ml-— " 3 , ^12:— " 3 1 " 



_^ ; 4 L' 

R 2 VI ' " 12 '~ 2R 2 vl "~ dT 2R 2 VV : 



_2V(T)uL 3 T 3 AV{T)ulLT dV(T)2u^LT 2 

13 : " rdI + VV + ~dT ~jv~' 

The four-dimensional mass matrix will get contributions from all the quadratic terms in 
5.2 Gauge fields 

Now consider the gauge fields and . We will use the symmetric and anti-symmetric 
combinations 

A ± {x^,u) = ^(A (1 \x",u)±A^ 2 \x^,u)), (5.4) 



10 In general there is also a pseudo-scalar fluctuation of the phase of the tachyon 8(x^, u), but we are working 
in unitary gauge where 8 = 0. 
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and consider only the fi and u components. The former will give rise to vector and axial- 
vector (A~) mesons, and the latter to scalar and pseudo-scalar mesons in four dimensions. 
Let us also further fix the gauge by setting 



At = 0. 

Expanding the action to quadratic order in the gauge fields then gives 



S[A + ,A~] = -M d 4 xduuV(T)Vv 



_ \p+ 2 , _\p- 2 _| ' 

2 1 ^' 2 1 ^ R 2 V 



+ 



4u2T 2 



l + u %T 2 L 2 



1 + 



2 



(5.5) 



Kl 2 + l^J 2 ) 



T 2 L 2 {L') 2 \ ,._, 2 4ufT 2 _ n2 2u^T 2 LL' 



i? 2 p 



A 2 P 



(5.6) 



5.2.1 The A + sector 

The action in the symmetric (vector) sector is given by 



S[A + ] = -N d 4 xduuV(T)Vv 



l + U 2T 2 L 2 



(5.7) 



This sector is similar to the gauge field in the single 8-brane case in ^ . We expand the gauge 
field A~^ in radial modes ip n (u) 



A+{x^u) = Y,4 {n) WMu), 



(5.8) 



that satsify the eigenvalue equation 

1 



V(T)uVV 

and the normalization condition 



a ( v(T f_±g2^, K 



R 2 VV 



M J duuV(T)Wip m ^ n = -5,, 
The four-dimensional action in this sector is then 
5[a+W] - 



(5.9) 



(5.10) 



where f^j^ = d^at^' - d v Op . 

The zero mode with = 0, is special. The eigenvalue equation ( ^9| ) gives 



dv- 



R 2 VV 



V{T(y)){v + v2T 2 L 2 ) 



(5.12) 
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In the UV asymptotic region this becomes 

V>o ~ / dv v~ 

J oo 



■21/2 



(5.13) 



so ipo(u) is UV-normalizable. On the other hand, in the IR asymptotic region we get 



ipo{u) ~ / dvexp- wT x < 

Juo (v ~ uoY 



v- 7 / A (v - n )~ 3 (p > 2) 

( v -7/4+ ^2 W -V4) ( v -u )- 3 (p = 2) (5.14) 

/^-l/4 (v _ no) 2 P -7_ (p<2) 



The exponential divergence of the integrand as v — > no implies that vpo{u) is non-normalizable 
in the IR. We therefore have to exclude this mode from the spectrum. The spectrum of vectors 
is therefore purely massive. Note that the exponential divergence is due to the tachyon. 



5.2.2 The A~ sector 

The action in the anti-symmetric sector is given by 

"1 



S[A~] = -AfJ d A xduuV(T)Vv 



where the coefficients are given by 



FZ\ 2 + B 1 \F-\' + B 2 \A-\' + B 3 (Azy + B A F-A-' i 



2"*' 



(5.15) 



_ l + u*T 2 L 2 

Bl ~ r 2 v ' 62 



Au^T 2 



1 + U 2T 2 L 2 



T 2 L 2 (L') 2 \ „ AuIt 2 „ 2uIt 2 LL> 

1 + ' 



AuiR 2 V 



,B 3 



R 2 V 



-,B 4 



R 2 V 



(5.16) 



We decompose the four-dimensional part of the gauge field A^ into a longitudinal component 
a\ and transverse components Aj^ (where d^A 1 -^ = 0), and expand all the components in 
radial modes 

A^,u) = 5>;("V)^(u), (5.17) 

n 

A\{x^u) = ^2d^ n \x^l(u), (5.18) 

n 

A-(x»,u)=^2J n \xnCn(u). (5.19) 

n 

These modes satisfy the eigenvalue equations 



d,, 



uV^VvUti'-^B^ 



uV(T)Vv (b 2 & - l -B A & - (m-) 2 ^) (5.20) 



Bs(n = M 2 n 



Bl(Cn-d') + 2 fi 4d 



(5.21) 
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and an additional equation relating the longitudinal and pseudo-scalar modes 

and the normalization conditions are given by 

M 



/ duuV(T)Vv 

J un 



/ duuV(T)Vv^ = -5 n 

JUQ Z 



+ 2 ^4 {(Cm - Cm)d + £m(Cn ~ £n)} 

The four dimensional action in this sector then becomes 



2 "ran • 



(5.22) 



(5.23) 



(5.24) 



S[a-to,uM] 



where = d^a v ^ - d^a^ 



n 
(n) 



(5.25) 



Using (|5.20|) , (|5.2l|) , ( |5.22D , (|5.23| ) and (|5.24D , we can express the masses m n and M n as 

(5.26) 
(5.27) 



(m-y = 2jV / duuV(T)VV B 1 (d u ^y+B 2 (^y-B^nd u ^ 



(M n ) 2 = 2J\f / duu^Ty^CCn) 1 



The spectrum in this sector consists of massive axial- vectors ^ and massive pseudo-scalars 
u>( n ). The pion is identified with the lowest pseudo-scalar mode Eq. ( |5.27 ) shows that 
the pion acquires mass due to the non-trivial tachyon background. Below we will estimate 
this mass for the case of a small quark mass. 

5.3 The pion mass 

We will begin by establishing that the pion is massless when the quark is massless. In this 
case the large u behavior of the tachyon and the 8-brane-anti-8-brane separation is given by 

T(u) « C T u- 2 e- RL °° u , L(u) m + C L u" 9 / 2 , (5.28) 

where Ct is related to the chiral condensate. For Mo = the solution to eq. (5.21) is 

Co(«) = 0. (5.29) 

The large u asymptotic form of ^|J( n ) can then be read off from the asymptotic behavior of 
eq. flfjp : 



« a + bu 3/2 , 



(5.30) 
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where a and b are constants which we will determine shortly. It can easily be checked that 
this solution is UV normalizable under the condition ( |5.24 ) , and therefore that it corresponds 
to a massless pion. 

Now turn on a small quark mass C' T . The condition for the validity of the corresponding 
tachyon solution 

T{u) u- 2 (C T e~ RLooU + C T e RL °° u ) , (5.31) 

is Ct ~ e^ RL °° u °° . We can therefore treat the quark mass as a perturbation of the massless 
solution we found above. To leading order in C T , and at large u, equation ( |5.21 ) gives 

Co ( u ) ~ (C' T e RL ™ u + C T e- RL ™ u y 2 , (5.32) 

where Mg = 0{C' T ). We have kept the C' T term from the asymptotic solution for the tachyon 
since it comes with a growing exponential. We can now use ( 5.27[) to express the pion mass 



Mo in terms of the parameters Ct and C' T . Since the integral in ( |5.27| ) is dominated by large 
u due to the presence of the tachyon potential, we get 

, r 2 9b 2 M f u °° du 

M n ~ — — / ^ , 5.33) 

° 8R 2 J U1 (C' T e RL °° u + C T e~ RL °° u ) 

where lower limit of the integral satisfies 1 <S u\ <C Uoo. Under the condition that C' T < 
e --Rioo«oo ) anc j assum i n g that Ct ~ 0(1), this gives 

Using ( p.13 ) and ( 3,1§| ), and inverting we get 



M o ~ l^ptf 171 ^) ■ ( 5 - 35 ) 

Let us now compute the constants a and b that appear in the longitudinal zero mode 
( |5~30| ) for zero quark mass. In principle, these constants are fixed by the boundary conditions 
and the normalization condition ( |5.24 ). However the latter requires more knowledge than 



we have about the form of the solution at finite u. Instead, what we will show is that these 
constants are related to the gauge theory parameter f n , the so-called pion decay constant. 
To see how, recall that in QCD the pion decay constant can be extracted from the two point 
function of the axial vector current in the massless quark limit 

fl = n A (0) = J d 4 x (J%(x) Ja^(O)) . (5.36) 

Using the usual holographic dictionary, we evaluate this by varying the action of the solution 
twice with respect to the boundary value of the dual field Aj^. Fourier-transforming in 
spacetime, the part of the action ( |5.15 ) that depends on this field becomes 



SIA 1 ] =N J ^duuV{T)^[BMtf + (P 2 +B2)\K\ 2 - B ± A t' A ^] » ( 5 - 37 ) 
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where the coefficients were defined in (|5.16| ) . The equation of motion reads 

d u uV^y/vUxA^-^BjA =uV(T)Vv^p 2 + B 2 )Aj I -^B 4 A^ . (5.38) 



Evaluating the action on the equation of motion then gives 
SiA 1 



M I S^du d 



(5.39) 



The IR boundary u = uq doesn't contribute since the tachyon potential V(T) goes to zero 
exponentially. Consider the zero mode A^{p^,u) = Co~( n ) a M ^ (P^)> an d impose the boundary 
condition Coi u oo) = 1- The pion decay constant is evaluated by varying with respect to a^/ ^ 
and imposing p 2 = 0: 



f 2 = NuV{T)Vv (b^u^u) - ~B4£<j-M 2 ) 



(5.40) 



To relate this to the longitudinal zero mode we note that the equation for £q (u) which follows 
from ( 5.3SD at p 2 = is precisely the same as the equation for £q(m) with (o(u) = 0. 

This implies that the two zero modes are proportional to each other. The proportionality 
factor can be determined from the normalization condition for £r,(it) (|5.24 ), which when 
combined with ( |5.22 ) gives 



i = N uV(T)Vv (b4( U )^'(u) - ~B4 



(u)' 



(5.41) 



Comparing with ( |5.40 ) we see that £q(u) = y/^f^^u). The boundary condition on the 



longitudinal zero-mode is therefore ^(""oo) = V(V2/t)> an d the condition ( 5.41 ) then fixes 
the other constant: 



1 V2f w R 2 _ 



V2f n SAf 



3/2 



(5.42) 



We finally get 



M, 



2 _ 4 m q(qq) 
f 2 

J 7T 



(5.43) 



which up to a factor of two reproduces the GOR relation. 
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6. Conclusions 



The Sakai-Sugimoto model provides a holographic description of a gauge theory that is close 
to large N c QCD with massless quarks. The most compelling feature of this model is that 
it exhibits spontaneous chiral symmetry breaking in a simple geometrical way. We have 
extended this model by adding the flavor bi-fundamental scalar field corresponding to the 
open strings between the 8-branes and anti-8-branes. This field is dual to the operator qq, 
and therefore describes both the quark mass deformation, as well as the chiral-symmetry- 
breaking quark condensate. 

Our analysis was carried out in the non-compact AHJK model, which is really dual to 
the NJL model. However, since the UV behavior in the compact and non-compact models 
is identical, our results for the quark mass, quark condensate, and GOR relation hold in the 
compact Sakai-Sugimoto model as well. The IR behavior of the solution will be quantitatively 
different, though qualitatively similar. 

In the non-compact model there exists also a parallel 8-brane-anti-8-brane configuration 
in which the chiral symmetry is unbroken. In this configuration the bi-fundamental field is 
tachyonic near the origin, so the configuration is locally unstable. The stable configuration is 
the connected, U-shaped configuration in which the chiral symmetry is broken. The same two 
configurations exist in the Sakai-Sugimoto (compact) model at high temperature. It would be 
interesting to study the effect of the bi-fundamental scalar in that case as well. In particular if 
the tachyonic mode becomes massless at some temperature, it might indicate a second-order 
(rather than first-order) phase transition to a chiral-symmetric phase. 

Far from the origin the bi-fundamental scalar field is massive, but it still has an impor- 
tant role in the holographic duality. The 8-brane-anti-8-brane theory with this field exhibits 
a U-shaped 8-brane solution with a non-trivial profile for the bi-fundamental field. At large 
distance (the UV of the gauge theory) the bi-fundamental field has a non-normalizable expo- 
nentially growing component, and a normalizable exponentially decreasing component. We 
showed that the former is related to the quark mass, and that the latter is related to the 
quark condensate. We have also found a numerical solution in the normalizable case, and 
compared it to the solution of the model without the bi-fundamental field. We showed that 
including this field lowers the free energy of the solution. 

Lastly, we began an analysis of the fluctuations of the 8-brane and anti-8-brane world- 
volume fields, including the bi-fundamental scalar, and the adjoint scalar and gauge fields. 
These correspond to the various mesons. In particular the lowest mode of the tt-component 
of the antisymmetric combination of the gauge fields describes the pseudo-scalar pion. We 
evaluated the mass of this mode in terms of the asymptotic behavior of the bi-fundamental 
field, and showed that it satisfies the GOR relation. In other words the pion mass is propor- 
tional to the product of the quark mass and the quark condensate. The rest of the meson 
spectrum should be analyzed numerically. 

There is an interesting question regarding the bi-fundamental field. In the action ( |2.lD 
the field T corresponds to a straight string along X4. However, as was shown in [11], the open 
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string between the 8-brane and anti-8-brane curves into a U-shape. It is not immediately 
clear to us how to incorporate this effect into the action. It would be interesting to study to 
what extent this affects the solution with T. 

Another interesting open question is the spectrum of fluctuations of the bi-fundamental 
scalar field T. In particular the spectrum of fluctuations around the parallel and U-shaped 
solutions with T = will determine their (in) stability. 
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